Drift wave turbulence in the presence of a dust density gradient 
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Abstract 

We present turbulent properties of electrostatic drift waves in a nonuniform collisional plasma 
composed of magnetised electrons and ions in the presence of immobile dust particles. For this 
purpose, we derive a pair of nonlinear quasi-two-dimensional equations exhibiting the coupling 
between the generalized ion vorticity and density fluctuations associated with collisional drift waves. 
The effect of a dust density gradient on the initial drift instability and fully developed turbulence 
is examined numerically. 
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I. INTRODUCTION 



Low-frequency (in comparison with the ion gyrofrequency) electrostatic drift waves in 
a nonuniform coUisionless magnetized electron-ion plasma are supported by inertialess 

n 

electrons and inertial ions [1]. In the pseudo-three dimensional drift wave electric field 
E = — V0, where is the wave potential, the electrons have a helical trajectory due to their 
= (c/-Bo)E X Bq and v^) = — (cA;BTe/ei?Qr2e)Bo x Vrig drift motions in a plane perpendic- 
ular to the external magnetic field Bq = zi?o, as well as their rapid motion along the external 
magnetic field direction, where z is the unit vector along the z axis in a Cartesian coordi- 
nate system, Bq the strength of the magnetic field, c the speed of light in vacuum, ks the 
Boltzmann constant, e the magnitude of the electron charge, T^. the electron temperature, 
and He the electron number density. 

Furthermore, the motion of the two-dimensional ions in a plane perpendicular to the axial 
direction is governed by the v^; and the polarization Vp = —{c/ BoUJci)dtV\(t) drifts, where 
we have denoted dt = dt + ^ e ■ ^ , = eBo/niiC is the ion gyrofrequency, and the ion 
mass. 

Hence, the electron and ion density fluctuations are different due to the differential elec- 
tron and ion motions in a magnetized plasma with an equilibrium density gradient. The 
resulting space charge separation, in turn, causes dispersive drift oscillations to propagate 
across the homogenous magnetic field and density gradient directions. 

In a coUisional electron-ion magnetoplasma, dispersive drift waves are excited by the 
combined action of the density inhomogeneity and electron-ion collisions. The latter produce 
a phase lag between the magnetic field-aligned electron velocity perturbation as well as the 
drift wave potential and the electron density perturbation. As a result, drift waves grow 
exponentially, extracting energy from the equilibrium density gradient. 

It is understood that nonthermal drift waves cause anomalous cross-field diffusion of the 



plasma particles 



-[3|. The plasma particle confinement is significantly improved if there 



are cylindrically symmetric zonal/sheared flows 



^-|6|] that are nonlinearly excited by finite 



amplitude dispersive drift waves in plasmas. Thus, zonal flows act like a transport barrier 

0. 
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Charged dust impurities are common in space and laboratory plasmas [8l4l0l|. The 
presence of charged dust grains in an electron-ion plasma modifies the equilibrium quasi- 
neutrality condition riio = ^eo + e-Z^d^dO; where rijo is the unperturbed number density of 
the particle species j [j equals i for ions, e for electrons, and d for charged dust grains), 

is the number of constant charges residing on dust, and e = 1 ( — 1) for negative (posi- 
tive) dust. Due to the modification of the quasi-neutrality condition, the difference between 
the divergence of the electron and ion fiuxes involving the V^; drift is finite. This leads 



to a fiute-like Shukla-Varma (SV) mode ll| in a nonuniform dusty magnetoplasma. The 



SV mode governs the dynamics of the ion vorticity, which evolves in the form of a dipolar 
vortex. 

In this paper, we present an investigation of the pseudo-three dimensional dissipative drift 
wave turbulence in a nonuniform coUisional dusty magnetoplasma composed of the electrons, 
ions and stationary charged dust grains. By using the two-fiuid model and the guiding 
center drift approximation (viz. \d/dt\ <C u^d), we derive a pair of nonlinear generalized 



Hasegawa-Wakatani (HW) equations [1^ governing the evolution of the potential and density 
fiuctuations associated with the low- frequency (in comparison with Ud), long wave length 
(in comparision with the thermal ion gyroradius) coUisional drift waves in our dusty plasma. 



Generalized HW type equations in the presence of dust have before been derived and 

13 1 and [l^. A linear analysis of dissipative drift modes in the 



linearly analyzed in refs. 



presence of a dust density gradient has been presented in Ref. [l5|, and further discussed in 



Ref. 



16| with respect to nonmodal growth, including dust charging effects. 



The work presented in refs. [I3l4l6l| (and others for similar topics) has shown that a 
dynamically variable dust charge can have a profound effect on linear drift wave stability. 
In a comprehensive dusty plasma turbulence model the charging process should therefore 
be taken into account. We have here refrained from including the charging process in order 
not to multiply the number of free parameters beyond the absolute necessary. 

The effect of dust charging dynamics on (nonlinear fully developed) drift wave turbulence 
has (to our knowledge) not yet been addressed in the literature. While we recognize that 
this would indeed be a worthwile effort, it is another kind of study itself 
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An important point why we have so far refrained from performing this study ourselves is 
the many uncertainties related to the dust grain charging model for a nonuniform, collisional 
magnetoplasma that have been used so far. In particular, the dust grain charging is a 
dyanmical process and requires a complete kinetic theory for fluctuating (in the electrostatic 
field of pseudo-three-dimensional drift waves) electron and ion currents that reach the surface 
of dust grains in our plasma. 

In the following, we study turbulence properties of the modified (by the presence of 
charged dust particles) drift waves in a nonuniform collisional magnetoplasma with static 
charged dust impurities. Thus, we are assuming that the timescales for the excitation of 
drift waves and their evolution are much shorter than the dust plasma and dust gyroperiods, 
and much longer than the electron-ion relaxation period. 

Furthermore, the assumption of constant dust charge would remain valid since the dust 
grain charging frequency is typically much larger than the modified drift wave frequency. 
Thus, our nonlinear simulation study of the dissipative drift wave turbulence is based on a 



slightly different set of dusty HW-SV equations compared to refs. |13l4l6|. excluding dust 
charge fiuctation effects. Our simulation results reveal that the equilibrium dust density gra- 
dient controls the formation of turbulent structures and the associated cross-field transport 
of the plasma particles in a nonuniform dusty magnetoplasma. 

The manuscript is organized in the following fashion. In Sec. II, we derive the governing 
nonlinear equations for the low-frequency (in comparision with the ion gyrofrequency) dissi- 
pative drift wave turbulence in a nonuniform dusty magnetoplasma composed of magnetized 
electrons and ions, as well as unmagnetized stationary charged dust grains. We use the two- 
fluid equations in the guiding center approximation to derive a pair of nonlinear equations 
that exhibits coupling between the plasma density perturbation and the drift wave potential. 
Section III presents a local linear dispersion relation and its analysis. Governing equations 
for nonlinearly interacting finite amplitude drift waves are numerically analyzed in Sec. IV. 
Summary and conclusions are contained in Sec. V. 
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II. NONLINEAR DRIFT WAVE EQUATIONS 



Let us consider the electrostatic drift waves with the electric field E = — V^^ — dzcf) in a 
plasma with a uniform magnetic field zi?o and the equilibrium density gradient d^njo. 

In a low-beta (/3 = SirneokBTf./ Bq ^ 1) plasma, the electron and ion fiuid velocities for 
\d/dt\ ^ Uci, i^e are, respectively, 

Ue — Z X V0 Z X Vne + ZMe^, (1) 

Bo eBoUe 

and 

with the magnetic field-aligned electron fiuid velocity perturbation is given by 

e d f kBTeUeA 
Uez = h (3) 



where rrie is the electron mass, z/g cjce = ^BQ/rriec) the electron-ion collision frequency, 
z/j the ion-neutral collision frequency, and riei = n^. — n^o ngo) is a small perturbation in 
the electron number density. 

The two-dimensional (2-D) ions are assumed to be cold. The assumption of 2-D ions 
ensures the decoupling of the dust ion-acoustic and drift waves 

Inserting eqs. [I]-|3]into the electron and ion continuity equations, we obtain the evolution 
equations for the electron and ion number density perturbations Uei rieo) and rin njo), 
respectively 

„^ -— zx Vn,o-V0 + :^— (0-^— ), (4) 

and 



dt Bq rrieVe dz"^ \ e Ueo 



- — z X Vnio ■ V0 - [-r^ + ^i] Vl0. (5) 



dt Bq BquJci \dt 

Now, subtracting eq. [5] from eq. H] we obtain the generalized ion vorticity equation 
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dt J dy z/e dz^ \ e n^o J 

1 /2 

where = (eZdc/i?o^jo)f^2.'^dO) ^lh = i^ce^d/d) is the lower-hybrid resonance frequency 
in our dusty plasma, the drift scale is ps = Cs/cod, and Cg = {akBTe/niiY^'^ is the modified 
ion-sound speed with a = nio/rieo > 1. 

We note that the term arises due to a non-zero value coming from the difference 
between the divergence of the electron and ion fluxes involving the (c/-Bo)E x z drift in our 
nonuniform dusty plasma In the absence of the magnetic field-aligned electron motion, Eq. [6] 



describes the dynamics of the S V mode ll| . In a quasi- neutral dusty plasma with immobile 
charged dust grains, eqs. [Hand [6] are closed with the help of rin = riei = ni, which is valid 
as long as the frequency (wavelength) of the drift waves is much larger than the plasma and 
dust gyrofrequencies (the electron Debye radius). 

It is appropriate to normalize the wave potential by ksTe/e, the density perturbation rii 
by Ueo, as well as the time and space variables by the ion gyroperiod uj~'^ and the effective ion 
sound gyroradius ps, respectively. Then, the governing nonlinear equations for the coUisional 
drift waves in our nonuniform dusty magnetoplasma are 

dn 86 f d \ ^ ^ 

where Q = n = ni/ueo, 7 = ^i/^^a, a = KdPs, b = KiPs, where Ki = 

— {l/neo)dniQ/dx > 0. We have Fourier decomposed the drift wave potential and elec- 
tron density perturbations along the z-axis, but kept their variations in the x — y plane, so 
that g = k'^Uce/i^eO: is the dissipative coupling coefficient. 



Equations and ([H]) are the generalization of the HW equations [12|]. The latter are 
recovered in the limit a = and 7 = 0. 
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III. LINEAR DISPERSION RELATION 



Linearization of eqs. ([7]) and ([8]) with an ansatz (j),n ^ exp[—iut + ■ x] obtains the 
dispersion relation 



u^k^ + ijA- 'jgk^ + iuoG - igB. (9) 

where A = aky, B = bky and G = gk"^ + '-yk'^ + g. Inserting oj = ujr + iuji and solving the 
imaginary part of the dispersion relation for uji ^ ujr results in the real frequency 

qB bky , . 

The growth rate is obtained from the real part as 

-1 



1 T VI + + 19)k'' + Aujr\ (11) 

with K = 2k/ G. Only the negative root is relevant for unstable solutions (w/ > 0). 
Instability thus occurs if 

{ujl + -ig)e + AujR>Q. (12) 

This is always achieved if Alor ~ AB ah > 0. The usual resistive drift wave instabilty 
(expressed by the first term) is thus enhanced for co-aligned gradients with a > if 6 > 0. 

The growth rate is reduced for counter-aligned gradients. Drift modes are damped (w/ < 
0) for a large counter-aligned dust density gradient, when \a\ > [k/{l + A;^)]^|6|. For k = 1 
this corresponds to a < —6/4. 

As drift wave turbulence is most active around k ^ 1, it can be expected that (for 6 = 1) 
the turbulence is strongly damped for a < —0.25, and enhanced approximately proportional 
to i/a for positive a 

In Fig. [Uthe growth rate ui{k, a) is plotted as a function of wave number and dust density 
gradient. 
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FIG. 1: Growth rate u}i{k,a) as a function of wave number k and dust density gradient a. Black 
marks the region of stabihty. 



IV. NONLINEAR NUMERICAL SOLUTIONS 

For numerical simulations, we write Eqs. ([7]) and (j8]) as 

dtn + [0, n] = -ady(j)--iVL + g {(p-n), (13) 
dtn+[(j),n] = -{b + a) dy(f) + g {(j) - n), (14) 

where the advection nonlinearity [A, B] = z x VA ■ VB is expressed in terms of the Pois- 
son brackets. The equations are numerically solved with a third order Karniadakis time 



stepping scheme in combination with the Arakawa method for the Poisson brackets |18N20|. 
Hyperviscuous operators z/^V^ with z/^ = —2 ■ 10^^ are added for numerical stability to 
the right-hand side of both equations ( |T3l) and ( |T4l) . acting on Q and n, respectively. The 
Poisson equation is solved spectrally. 

The equations are for the turbulence computations discretized on a doubly periodic 512 
X 512 grid with box dimension = Ly = QApg. The initial drift wave development are 
computed on a 128 x 256 grid with box dimension = 64ps and Ly = 128ps Nominal 
plasma parameters are g = 0.5 and 6 = 1. 
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V. VORTEX DEVELOPMENT 



First, we study the initial evolution of a drift wave growing out of a density perturbation 
in the presence of a dust density gradient a. 

The ion viscosity is set to 7 = 0.1. A Gaussian density perturbation is initialized with 
amplitude n = 1, which generates a potential perturbation of the same initally Gaussian 
form This leads to the formation of an E x Bq drift vortex, where the plasma is convected 
around the perturbation. 

For a = the usual drift wave dynamics is encountered, which is shown in Fig. [2^: the 
singular density structure {t = 0, left frame) is propagating in the electron diamagnetic drift 
direction (upwards, in the positive y direction) and generating drift wave structures in its 
wake, which are shown for times t = 25 (middle frame) and t = 50 (right frame). 

At a later stage, the drift wave nonlinearly develops via secondary instabilities into a 
saturated turbulent state, whose properties are discussed in the next section. 

In Fig. the evolution for a = —0.5 is shown: the vortex is elongated in the y direction, 
while the structure is propagating but damped. The mode is still damped when the viscosity 
7 is set to zero and the initial amplitude to no = 10. 

Figure [2]d shows the case of a = +0.5: The structure is much slower propagating in the 
y direction, but is forming a nearly standing drift wave structure, which is elongated in the 
X direction, and growing fast in amplitude until nonlinearly secondary modes form by the 
Kelvin-Helmholtz break-up of the streamers. 

The linear and nonlinear development of the drift wave perturbations thus strongly de- 
pends on the sign of the dust density gradient. Co- alignment of the plasma density gradient 
with the dust density gradient is enhancing the mode growth, while counter-alignment is 
strongly damping. 
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FIG. 2: Development of a drift vortex with an initial Gaussian density perturbation n(x, y) in the 
presence of a dust density gradient a for t = (left), t = 25 (middle) and t = 50 (right): (a) 
a = 0, (b) a = —0.5, (c) a = +0.5. Note: the color scale (red for positive, blue for negative density 
perturbations in drift scaled units) for n(x,y) in case (c) is 30 times larger than in (a) and (b). 
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FIG. 3: The potential energy U (hold), the kinetic energy K (thin) and the generalized enstrophy 
W (dashed) as a function of the dust density gradient a. 

VI. TURBULENT STATE 

The computations are started from a random quasi-turbulent spectral bath and run into 
saturation. The spectral bath is generated initialising the density field in Fourier space by 
(within some range) random amplitudes that approximately follow a power law spectrum. 
In this way all k modes are excited from start, and a fully developed turbulent state is 
reached rather rapidly in the simulation. 

Statistical averages are taken in the interval between 200 < t < 1000, which for all 
parameters is well in the saturated regime. 

The turbulent activity can be characterised by the potential energy f/ = (1/2) J dVii"^, the 
kinetic energy K = (1/2) J dV{'V(j)y and the generalized enstrophy W = (1/2) J dV{n — 
Q)"^. The mean of these energetic quantities together with the standard deviation of the 
fiuctuations in the saturated state are displayed in Fig. [3] for various values of the dust 
gradient parameter a. 

The turbulent energies closely follow the predictions from the linear analysis: negative 
values of a (counter-aligned gradients) have a strong damping infiuence on the turbulent 
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FIG. 4: The wave number spectra of density Euctuations in the x direction, multiplied hy a factor 
[kxPsY'- a = (hold line), a = \, a = —0.2 (dashed line). 

fluctuations. Co-aligned gradients are enhancing the drive of the drift wave turbulence. 

The wave number spectra \n{pskx)\{pskx)^ are shown in Fig. HI The standard HW case 
(a = 0, bold line) is arbitrarily scaled to 1. For a co-aligned dust density gradient (a = 1, 
thin line) modes around 0.5 < psk^ < 1 in the order of a few drift scales ps are in relation 
more strongly pronounced. For counter- aligned gradients (a = —0.2, dashed line) the smaller 
scales are most strongly damped. 

This change in structural scales is also visible in 2-D x — y plots of density perturbations 
(left column) and vorticity (right column) in Fig. |5l The top row shows the standard case 
for a = of the typical drift wave turbulence. In the middle row, a strong co-aligned 
dust density gradient a = 1 leads to a stronger perturbation of structures at smaller scales, 
visible by a more frayed out density field (contour line drawn for n = 0) and more strongly 
pronounced small-scale vorticity structures. The bottom row shows the case of a = —0.2 with 
counter-aligned gradients, where quasi-linear drift wave structures appear, which extend in 
the x-direction and propagate in the y-direction. 
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FIG. 5: Density n{x,y) (left) and vorticity Q{x,y) (right) for (a) a = 0, (b) a = 1, (c) a = —0.2. 
Only a quarter of the actual computational domain is shown here. 
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VII. SUMMARY AND DISCUSSIONS 



Summarizing, we have investigated the properties of nonhnearly interacting finite amph- 
tude electrostatic drift waves in a nonuniform coUisional dusty magnetised plasma composed 
of 3-D magnetized electrons, 2-D magnetized ions, and immobile charged dust impurities. 
The dynamics of the dissipative drift wave turbulence in our model of a dusty magnetised 
plasma is governed by generalized Hasegawa-Wakatani equations In the latter, one encoun- 
ters nonlinear coupling between the quasi-neutral electron density fluctuation and the ion 
vorticity that are driven by the combined action of the dust density gradient and a dissipa- 
tive electron current arising from the magnetic field-aligned electron motion reinforced by 
the parallel electric force and the parallel variation of the electron pressure perturbation for 
the constant electron temperature. 

Numerical simulations of the governing nonlinear equations reveal features of the fully 
developed drift wave turbulence The presence of a co-aligned dust density gradient is found 
to enhance the drift wave turbulence and the cross-field plasma particle transport. On the 
other hand, counter-aligned gradients lead to a damping of the drift fluctuations. Thus, 
the dust density inhomogcneity plays a decisive role in the formation of drift wave vortex 
structures in nonuniform dusty magnetised plasmas 

In conclusion, we stress that the present results should be helpful in understanding the 
features of fully developed low-frequency drift wave turbulence that might emerge from the 
forthcoming laboratory dusty plasma experiments in an external magnetic field, and from 
nonuniform magnetized space dusty plasmas. 

Acknowledgements 

The work of AK was supported by the Austrian Science Fund (FWF) Y398 and by 
a junior research grant from University of Innsbruck. The research of PKS at UCSD was 
supported by NSF PHY0903808. The authors thank M. Rosenberg (Department of Electrical 
Engineering, University of California San Diego) for valuable discussions. 



14 



[1] B. B. Kadomtsev, Plasma Turbulence (Academic, New York, 1965). 
[2] W. Horton, Rev. Mod. Phys. 71, 735 (1999). 

[3] J. Weiland, Collective Modes in Inhomogneous Plasma (Institute of Physics, Bristol, 2000). 
[4] P. K. Shukla, M. Y. Yu, H. U. Rahman et al, Phys. Rev. A 23, 321 (1981); Phys. Rep. 105, 
229 (1984). 

[5] P. K. Shukla and L. Stenflo, Europhys. J. D 20, 103 (2002). 

[6] P. H. Diamond, S. I. Itoh, K. Itoh, and T. S. Hahm, Plasma Phys. Control. Fusion 47 R35 
(2005). 

[7] G. R. Tynan, A. Fujisawa, and C. McKee, Plasma Phhys. Control. Fuson 51, 11301 (2009). 
[8] P. K. Shukla and A. A. Mamun, Introduction to Dusty Plasma Physics (Institute of Physics, 
Bristol, 2002). 

[9] P. K. Shukla and B. Eliasson, Rev. Mod. Phys. 81, 25 (2009). 
[10] S. I. Krasheninnikov, R. D. Smirnov, and D. L. Rudakov, Plasma Phys. Control. Fusion 53, 
083001 (2011). 

[11] P. K. Shukla and R. K. Varma, Phys. Fluids B 5, 236 (1993). 
[12] A. Hasegawa and M. Wakatani, Phys. Rev. Lett. 50, 682 (1983). 

[13] S. Benkadda, P. Gabbai, V. N. Tsytovich, and A. Varga, Phys. Rev. E 53, 2717 (1996). 

[14] S. V. Vladimirov and V. N. Tsytovich, Phys. Rev. E 58, 2415 (1998). 

[15] S. Benkadda and V. N. Tsytovich, Plasma Physics Reports 28, 395 (2002). 

[16] P. Manz and F. Greiner, Phys. Plasmas 17, 063703 (2010). 

[17] P. K. Shukla and V. P. Silin, Phys. Scr. 45, 508 (1992). 

[18] G.E. Karniadakis et al, J. Comput. Phys. 97, 414 (1991). 

[19] A. Arakawa, J. Comput. Phys. 1, 119 (1966). 

[20] V. Naulin and A. Nielsen, SIAM J. Sci Comput. 25, 104 (2003). 



15 



